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ABSTRACT 

Since the launch of the Fermi Gamma — ray Space Telescope on June 11, 2008, 55 gamma-ray 
bursts (GRBs) have been observed at coordinates that fall within 66° of the Fermi Large Area 
Telescope (LAT) boresight with precise localizations provided by the NASA Swift mission or other 
satellites. Imposing selection cuts to exclude low Galactic latitudes and high zenith angles reduces 
the sample size to 41. Using matched filter techniques, the Fermi /LAT photon data for these fields 
have been examined for evidence of bursts that have so far evaded detection at energies above 100 
MeV. Following comparisons with similar random background fields, two events, GRB 080905A and 
GRB 091208B, stand out as excellent candidates for such an identification. After excluding the six 
bright bursts previously reported by the LAT team, the remaining 35 events exhibit an excess of LAT 
"diffuse" photons with a statistical significance greater than 2 a, independent of the matched filter 
analysis. After accounting for the total number of photons in the well-localized fields and including 
estimates of detection efficiency, one concludes that somewhere in the range of 11% to 19% of all 
GRBs within the LAT field of view illuminate the detector with two or more energetic photons. These 
are the most stringent estimates of the high energy photon content of GRBs to date. The two new 
events associated with high energy photon emission have similar ratios of high to low energy fluences 
as observed previously. This separates them from bursts with similar low energy fluences by a factor 
of ten, suggesting a distinct class of events rather than a smooth continuum. 
Subject headings: gamma rays: bursts 



1. INTRODUCTION 

One of the more surprising results of the Compton 
Gamma-Ray Observatory was the EGRET discovery of 
an 18 GeV photon associated with GRB 940217 (Hur- 
ley et al. 1994). This was not an entirely unique event; 
about a half-dozen bursts were seen over the course of 
the mission with photons above 100 MeV (Catelli et al. 
1998, Dingus 2003). Since the GRB spectral energy dis- 
tribution at lower energies has been well characterized by 
a modified power law with peak fluxes at energies of the 
order of 200 KeV, the existence of photons at energies 10 4 
times higher puts a significant constraint on any viable 
model of the GRB phenomenon (e.g. Band et al. 2009). 
The obvious question is whether high energy photons are 
associated with all GRBs or only with a small sub-class. 
We are now in a unique period when two space missions 
are providing GRB coordinates with arc-second accuracy 
(Swift; Gehrels et al. 2004) and the ability to measure 
photon energies to 300 GeV (Fermi /LAT; Atwood et al. 
2009). Once these two missions cease functioning, the 
resources to investigate GRBs over a broad energy range 
are not likely to be available for decades. Thus, it is a 
matter of some urgency to extract the maximum amount 
of information from the assets now at hand. Prior to 
launch of the Fermi Gamma — ray Space Telescope, 
it was possible to speculate that the LAT instrument 
would detect more than 200 GRB events per year (Din- 
gus 2003). In the two-year period since the launch of 
the Fermi Gamma-ray Space Telescope, the Gamma- 
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ray Burst Monitor (GBM; Meegan et al. 2009) has 
reported approximately 475 GRBs, ie. a rate of about 
250 per year. Over essentially the same period, only 17 
bursts 3 have been identified by the Fermi /LAT. Given 
that roughly 50% of all GBM triggers point to directions 
within 66° of the LAT boresight, the fraction of bursts 
with photons above 100 MeV is apparently not much 
larger than 7%. Thus, it is an important task to apply 
the best possible techniques to explore whether this frac- 
tion would be significantly bigger if the detection thresh- 
old could be lowered. Such studies might significantly 
inform the design of future missions while extracting ad- 
ditional information from rather expensive facilities. 

From launch through March 16, 2010, 55 GRBs have 
been detected (mainly by Swift, as well as INTEGRAL; 
Winkler et al. 2003, and AGILE; Giuliani et al. 2008) 
at celestial coordinates that were simultaneously viewed 
by the Fermi /LAT. The angular resolutions of the Swift 
BAT and XRT are arc-minutes and arc-seconds respec- 
tively. Compared with the mean point spread width for 
the Fermi /LAT of a few degrees, these GRB coordinates 
are known with zero error. This sample is an excellent 
target for statistical methods that can take advantage 
of the precisely determined source direction. The basic 
technique that is employed here is the matched filter, 
most familiar to those detecting signals in the time do- 
main. The underlying assumption of this method is that 
the characteristics of both the signal and background are 
a priori known functions of one or more variables. Since 
the matched filter maximizes the signal-to-noise ratio, 
moderate departures from optimality degrade the filter 
performance relatively slowly, making this a valuable tool 

3 http: / / fermi . gsfc. nasa.gov / ssc / observations / types /grbs/ grb_table / 
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for investigating the possible existence of faint signals. 
We present the sample and data processing in the next 
section, and details of the signal detection technique in 
Section 3. Results are given in Section 4 followed by a 
brief discussion in Section 5. 

2. SAMPLE SELECTION AND DATA PROCESSING 

A selection of 148 precisely located GRBs (<1.8') was 
obtained from the Swift Web page 4 devoted to cata- 
loging all bursts. The data correspond to the period 
from launch until the end of March 2010. Most of these 
events were observed by the Swift/XRT, augmented by a 
few from INTEGRAL and AGILE. The resulting list was 
cross-correlated with Fermi spacecraft attitude data ob- 
tained from the Fermi Science Support Center 5 (FSSC) 
to identify only those events within 66° of the LAT bore- 
sight. 

Having identified the principal focus of this investiga- 
tion, two sets of LAT photon data were obtained from 
the FSSC corresponding to the 55 precisely located GRB 
fields and to 464 additional fields taken at random on 
the sky to study the background behavior. Through- 
out the following, they will be identified as "GRB" or 
"random" data. The FSSC data selection criteria cor- 
responded to any of the three event classes, "transient" , 
"source" or "diffuse" , photon energies above 100 MeV 
but below 300 GeV, a zenith angle less than 105° and 
time intervals when the DATA.QUAL parameter was 1 
and the South Atlantic Anomaly parameter, IN_SAA, 
was 0. The 105° zenith angle limit might be considered 
a bit risky but with an orbiting altitude of 565 km, the 
Earth's limb appears at a zenith angle of 112.83° which 
is easily resolvable. For the random LAT fields, fictitious 
GRB directions were chosen randomly over the sky but 
constrained by the less than 66° boresight requirement. 
Further cuts were imposed of less than 10.5° for the an- 
gle between the photon and GRB directions and a time 
window spanning from to — 100 s to to + 150 s. With 
LAT point spread width maxima of the order of 4°, the 
10.5° cut was imposed to allow up to 2.6 a errors without 
worrying about the effects of significant variations in the 
background flux over the effective solid angle of 0.105 sr. 
All data fields were additionally constrained to lie above 
or below the Galactic plane by more than 10°. Such ad- 
ditional criteria reduced the number of GRB fields from 
55 to 41. Six of these have been previously associated 
with significant fluxes of energetic photons (> 100 MeV) 
and thus provided guidance about the characteristics of 
the signals we were seeking. 

3. SIGNAL DETECTION TECHNIQUE 

The central theme of this paper is the use of a matched 
filter to enhance the signal-to-noise ratio for the detec- 
tion of GRB photons in a diffuse background. There 
are four independent variables that are involved: energy 
(E), photon angle with respect to the GRB direction 
{&), photon time of arrival (t) and photon event class 
(c G {1,2,3}). The most obvious consideration is that 
the photon direction should be correlated with the GRB 
coordinates within an angle related to the LAT point 
spread function for that specific photon energy. Pho- 



tons that fall far outside this criterion are almost cer- 
tainly background. The energy criterion is less intuitive. 
The energy spectra for both the GRB signal and the 
background are approximately represented by decreasing 
power laws with exponents Tqrb and T^ack respectively. 
Unfortunately, the background has a somewhat harder 
spectrum with T^ack ~ 1.5 while Tqrb — 2.2. Thus, 
higher energy photons are more likely to represent back- 
ground than signal. Photon directional accuracy is also 
critical so photons with small uncertainty (ie. high en- 
ergy) should be weighted more heavily. The time filter 
is the most problematic because GRB light curves lack a 
uniform shape. Particularly for the most energetic burst 
photons, the time delay relative to the low energy trigger 
can be substantial, up to thousands of seconds (Hurley 
et al. 1994). However uniquely identifying such late 
quanta with a GRB can be problematic, especially in 
view of the lack of a distinguishing spectral characteris- 
tic. One feature that is generally observed is a prompt 
flux within a few seconds of the low energy trigger. This 
can be seen in Figure 1 of Ghisellini et al. (2010). Most 
of the depicted LAT light curves exhibit a significant flux 
within 10 s of the burst. There are a few exceptions such 
as GRB 090328 so it must be recognized that not every 
fish can be caught with the same net. This led us to 
adopt a standard GRB light curve described by a bro- 
ken power-law with a break at t\, — 5 s. Prior to the 
flux rises as (t/tb) a ° after which it falls as (t/tb)~ a - The 
value for ao was set rather arbitrarily at 1 /2 while a was 
picked to be 7/5 in concert with the general behavior of a 
large number of GRB afterglows. We emphasize that all 
the filter parameters were determined prior to applying 
them to either the GRB or random data sets. 

When we began this analysis, we were unaware of 
the nuanced distinctions between the three LAT photon 
event classes, "transient", "source" and "diffuse", iden- 
tified as 1, 2 and 3. Since we were looking for a tran- 
sient point source, it seemed reasonable to include the 
"transient" quanta on equal footing with the other two. 
Only after going through the complete development of 
the matched filter search did we discover that the signal 
content of the "transient" photons is extremely small, 
~ 3.5%. This was determined by looking at the event 
class photon counts for the most intense GRBs, 080916C, 
090510, 090902B and 091003. The ratio, r GRB , of "tran- 
sient" and "source" GRB photons to "diffuse" is 0.34 and 
0.44 respectively, roughly independent of the particular 
burst. On the other hand, the same ratios for the se- 
lected random LAT fields, r^ack, are 9.65 and 0.85. Con- 
sequently, we realized that we needed to include a weight- 
ing factor, w c , based on event class that reflected this 
wide dispersion of information content. Serendiptously, 
the simple existence of a strong correlation of GRBs with 
"diffuse" class photons provided an independent verifica- 
tion of the validity of the matched filter technique. 

From these considerations, the mathematical form of 
the matched filter weight functions is given by: 

1 E 
Aiza 2 {E) E th 



4 http://swift.gsfc.nasa.gov/docs/swift/archive/grb_tablc.html/ 

5 http: / /fcrmi. gsfc.nasa.gov/ssc/ 
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rc-(±) a ° 0<t<t b 

Wt = \ c -(r b )~ a h<t<t c (3) 

I c = 2.056382; t c = ft b , 

-c= (Si)/(^ffi) 1<^<3 (4) 
w(E,6,t,c) = w E -wo -WfW c . (5) 

The point spread function for the LAT has been specified 
by the FSSC team in terms of a half-angle that includes 
68% of a random sample. We have fit this to a power 
law in energy (Akerlof & Yuan 2007) to obtain separate 
parameter sets corresponding to conversion in either half 
of the LAT tracker using the form: 

*66%{E) = (-£-)-**(> + vi, (6) 

For this analysis, we have assumed that the LAT PSF is 
Gaussian as described above for we with a characteristic 
width determined by the containment condition: 

<jpsf(E) = ^a 68% (E). (7) 

Although the actual LAT angular error distribution func- 
tion exhibits non-Gaussian tails, this approximation is 
unlikely to affect subsequent conclusions. As indicated 
above, the weight value for each photon is computed as 
a product of the four elements, we, we, wt and w c . The 
time filtering algorithm substantially ignores all photons 
that occur more than a few tens of seconds after the GRB 
trigger. For the parameters selected, this can best be de- 
scribed by an effective T go window time for the signal of 
16.864 s starting 1.624 s after the burst. The sum over 
all photons, Y] Wj, is the statistic used to identify events 
with significant fluxes of high energy photons. In as- 
sessing this technique, the mean contribution from back- 
ground photons turns out to be negligibly small. The 
identification problem is driven entirely by infrequent but 
large fluctuations of ^2 Wi contributed by background 
quanta. As an additional safeguard, a weight-sharing 
measure was introduced to ensure that the entire event 
weight was not generated by a single extreme photon. 
This parameter is defined by: 

= 3(w -wi -W2) 1/3 , g . 

W + Wl + W 2 

In this formula, wq, w\ and W2 are the maximum three 
photon weights for a particular event. The combined 
statistic, (J2 w i, incorporates both the sum of pho- 
ton weights and the additional requirement that no sin- 
gle photon contributes an overwhelming fraction of the 
value. The downside is that at least three photons must 
be associated with an event, an issue that will arise later. 

The significances of the Wi and statistics 
were determined by Monte Carlo simulation of the back- 
ground. For each field of the GRB data set, 2000 events 
were generated using the measured background spectral 
distribution to determine energy and uniform distribu- 
tions over solid angle and time to determine 9 and t. 
The number of background photons for each field was 
set by the measured background rate. By looking at the 
weight statistics for simulated events with and without 
injected GRBs, we gained reasonable confidence that the 



matched filter algorithm would correctly identify a sig- 
nificant fraction of real bursts. This work also demon- 
strated that the modified sharing-compensated statistic, 
(^2wi, was the stronger measure. 

As described earlier, 41 GRB fields were selected for 
examination. These are listed in Table 2. Six of these 
events have been previously identified with >100 MeV 
LAT-detected photons. These are GRB 080916C (Abdo 
ct al. 2009a), GRB 090323 (Ghiscllini et al. 2010), GRB 
090328A (Ghisellini et al. 2010), GRB 090510 (Abdo 
et al. 2009b), GRB 090902B (Abdo et al. 2009c) and 
GRB 091003 (Ghisellini et al. 2010). 

The 35 GRB fields with no previous claims for LAT 
detections were the target of this investigation. We rec- 
ognized that the most convincing argument for a true 
LAT identification should rely on the statistical distri- 
butions for the matched filter weights in LAT fields with 
similar characteristics. To increase that number as much 
as possible, the LAT fields of view were segmented into 24 
circular tiles embedded on a spherical surface as shown in 
Figure 1. Each tile subtends a cone with a half- angle of 
10.5°. This tiling scheme was applied to both the GRB 
and random field data sets to realize 697 and 8802 in- 
dependent directions in space satisfying all the criteria 
described previously. Taking advantage of the fact that 
each field observation was blocked into a 250-s segment, 
the number of independent observations was multiplied 
by five by regarding each 50-s time slice as a separate 
sample. Thus, there are 3485 background measurements 
taken from LAT observations obtained simultaneously 
with the candidate GRB fields and an additional 44010 
samples taken under similar but not identical conditions. 
One particular concern for an analysis of this type is that 
false positives will selectively occur as the sample photon 
rate rises substantially above the mean. Evidence that 
this is not the case here is shown in Figure 2 which plots 
the cumulative distribution of the total number of pho- 
tons within the LAT FoV over a 250-s interval. These 
rates explicitly exclude contamination from photons be- 
yond the 105° zenith angle cut. As shown in the plot, 
the two low fluence GRB events reported here are not as- 
sociated specifically with fields with high ambient back- 
ground rates. The similarity of the distributions for GRB 
and random fields also shows that the GRB data are not 
pathological as far as rates are concerned. 

4. RESULTS 

The statistical weighting scheme with the additional 
fractional sharing requirement clearly identified four of 
the six events previously associated with high energy 
gamma-ray emission. The exceptions were GRB 090323 
and GRB 0903028A for which photon counts were very 
sparse during the first 50 seconds following the burst trig- 
ger. Of the two, 090323 was marginally detected by the 
matched filter technique but GRB 0903028A provided 
no sign of a signal whatsoever. This gives some order 
of magnitude indication of the likely efficiency of the 
method. Three of the six were extremely intense with 100 
or more energetic photons (GRB 080916C, GRB 090510 
and GRB 090902B). As will be discussed later, the ex- 
istence of such a substantial bright fraction is a rather 
surprising situation. 

Most significantly, our analysis identified two addi- 
tional bursts, GRB 080905A and GRB 091208B, as prob- 
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Fig. 1.— A 10.5° tiling of the LAT field of view. 24 disks 
are closely confined within the 68° radius of the dotted circle that 
characterizes the acceptance of the LAT instrument. 




400 600 800 1000 1200 1400 1600 
number of photons 

Fig. 2. — Cumulative distributions of LAT photon rates over 
250-s intervals for the GRB (blue) and random background (red) 
fields. These rates reflect the entire LAT FoV except for photons 
that lie outside the 105° zenith angle cut. The rates corresponding 
to the two claimed GRB detections are indicated by arrows. 

ably associated with energetic photon detections in the 
LAT. The best estimate for the probability of these 
occuring by chance alone was obtained by performing 
identical searches on random LAT fields with the same 
criteria. Thus, 1 out of 44010 random fields gener- 
ated matched filter weights exceeding the value for GRB 
080905 A and 11 similar fields generated weights exceed- 
ing the value for GRB 091208B. Multiplying by a trials 
factor of 35 for the number of localized fields considered, 
the associated probabilities are 8.0 x 10 -4 and 8.7 x 10~ 3 . 
To check that these correlations were simply not due to 
preferentially higher background rates for the GRB ex- 



posures, we also performed similar calculations for the 
LAT data confined to ~ 20 uncorrelated directions and 
five independent time intervals from the same GRB data 
sets. All of these results are plotted as cumulative dis- 
tributions in Figure 3. The statistical similarity of the 
GRB off-axis and random field data demonstrates that 
the GRB data set is not correlated with anomalous en- 
vironmental conditions such as higher cosmic ray back- 
ground rates. A list of photons associated with these two 
events is provided in Table 1. 

Two independent observations at significance levels of 
0.1% and 1% already make this a reasonably credible re- 
sult. However, as we discovered the importance of LAT 
"diffuse" photons for identifying point sources such as 
GRBs, we realized that independent of the matched fil- 
ter weight analysis, a simple count of the number of such 
quanta was statistically quite significant. Thus for the 
35 well-localized GRB fields that don't include the previ- 
ously reported LAT-detected GRBs, the total number of 
"diffuse" photons is 25 for a rate of 0.714±0.143. For the 
GRB off-axis fields and the random fields, the compara- 
ble numbers are 1384 for 3485 fields and 17414 for 44010 
to yield per field rates of 0.397 ±0.011 and 0.396 ±0.003. 
The random probability of such differences is less than 
2.7%. As can be seen from the lists of photon counts 
in Tables 1 and 2, the two new GRB identifications are 
strongly correlated with these "diffuse" photons as well. 
Thus, the combined probablility that these claims might 
be randomly generated must be less than 1 x 10~ 5 . 

GRB 080905A is a short burst with optical counter- 
parts and a redshift of 0.128 (Rowlinson et al. 2010) 
while GRB 091208B is a long burst also with optical 
counterparts (Pagani et al. 2010) and a redshift of 1.063 
(Wiersema et al. 2010; Perley et al. 2010). As listed 
in Table 1, the estimated number of >100 MeV photons 
for each event is about 3. This is substantially lower 
than any previous LAT GRB identifications. Estimating 
the corresponding fluence of >100 MeV photons raised 
some interesting statistical issues that are perhaps not 
well appreciated. To first order, one would anticipate 
simply adding up all the energies of the detected quanta 
and dividing by the effective detector area. Although this 
produces a number, it bears only a casual relation to the 
intrinsic properties of the GRB in question. The reason 
is that the spectral number distribution, dN/dE follows 
a power law with exponent between -3 and -2. Thus the 
mathematical mean is defined but the variance is infi- 
nite. This leads to the conclusion that sums of energies 
drawn from the identical parent population will vary by 
factors of the order of unity, independent of the number 
of photons in the sample. In our case, the statistical un- 
certainty is already substantial since only 3 photons are 
detected for each event. Thus, we opted for a solution 
that invoked the least number of new assumptions: we 
adopted the LAT fluence values for bright bursts com- 
puted by Ghisellini et al. (2010) and listed in their Ta- 
ble 1. We next assumed that the spectral distribution 
for all bursts was an approximately universal function of 
energy so that fluences would scale with photon num- 
ber. This led to estimates of the LAT fluence value of 
3 x 10~ 6 ergs /cm 2 for both events. Unfortunately, the 
proportionality of fluence to photon counts is not very ev- 
ident. As can be seen in Table 2 for the two bursts, GRB 
090510 and GRB 090902B, for essentially the same num- 
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Fig. 3. — Complements of the cumulative distributions for C"l2 w i 
for 35 well-localized GRB fields (blue), 3485 random fields obtained 
nearly simultaneously with the GRB data (red) and 44010 random 
fields obtained at random times (green). 



TABLE 1 

List of high energy photons for GRB 
080905A & GRB 091208B 



GRB 080905A f 


= 0.2118 




= 69.66 


i t (s) 


(°) 


E (MeV) 


c 


Wi 


1 10.019 


0.21 


474 


3 


267.07 


2 5.661 


2.27 


181 


3 


61.06 


3 6.727 


10.41 


103 


2 


0.77 


GRB 091208B C 


= 0.8442 




= 17.23 


i t (s) 


(°) 


E (MeV) 


c 


Wi 


1 17.932 


2.06 


131 


3 


11.18 


2 20.985 


2.72 


125 


2 


3.94 


3 4.968 


1.66 


544 


1 


3.02 


4 6.364 


1.16 


1179 


1 


2.07 


5 19.735 


7.17 


197 


3 


0.20 



ber of photons, the fluences differ by a factor of 16. This 
problem may require a better realization of how fluences 
should be estimated. For future purposes when dealing 
with events with more numerous quanta, there is a ro- 
bust measure of fluence: multiply the number of photons 
by the median energy. In the continuum limit, a power- 
law distribution generates the following relation between 
Emean and E me di a n over the interval from threshold to 
infinity: 

Emean — p ^ ^T — 1 Emedian (9) 

Thus a robust proxy for gamma-ray burst fluence can be 
obtained by multiplying this estimate of E mean by the 
number of detected photons. 

The fact that >100 MeV photons have been observed 
with both long and short bursts suggests that the pro- 
duction mechanism is generic to relativistic jets, indepen- 



dent of the detailed structure of the progenitor system. 
So far, most of the high energy photon observations have 
been associated with large fluences at the lower energies 
measured by the Swift/BAT and the F 'er mi /GBM de- 
tectors. To explore such fluence correlations, Figure 4 
shows a scatterplot of LAT fluences vs. GBM fluences 
for the six bright bursts previously identified by the LAT 
collaboration, the two bursts uncovered by the statistical 
analysis described in this paper and the co-added sum of 
nine bursts with no strong perceptible radiation above 
100 MeV. For these purposes, GBM fluences were de- 
fined by photons in the range of 8 to 1000 KeV within 
the T90 time interval. The LAT fluences were obtained 
from Ghiselini et al. (2010) and the scaling method de- 
scribed previously. 

The salient feature of Figure 4 is the overall linear cor- 
relation of the LAT and GBM fluences with a ratio of 
unity between the two. The short burst identified by 
our analysis, GRB 080905A, lies above the trend line, a 
feature previously noted for this class of events (see dis- 
cussion of GRB 081024B in Abdo et al. 2010 and GRB 
090510 in Ackermannet al. 2010). GRB 091208B is com- 
pletely consistent with the higher fluence events. An in- 
teresting aspect of this graph is the estimate for the total 
co-added LAT fluence for nine bursts with GBM fluences 
roughly equal or greater to GRB 091208B but which lie 
a factor of ten lower but with rather large errors. If one 
believes that all GRBs have a similar generic behavior, 
it is possible that a large fraction of these events have 
delayed high energy emission that evades our estimates. 
It is more likely that there is a real population difference 
between events like 080916C, 090926A and 091208B and 
the majority of all other GRBs. Examining this ques- 
tion with larger data sets may shed some light on the 
mechanisms that drive high energy emission processes. 

Although the statistics are admittedly crude, we can 
use this analysis to set some limits on the fraction of high 
energy photon GRBs at the present levels of sensitivity. 
Starting with a sample of 41 well-localized events, wc 
have identiifed two with high statistical significance. The 
remaining number of LAT "source" and "diffuse" pho- 
tons above background is no more than approximately 
three. The efficiency of this technique is about 50% for 
a total so as many as 16 GRBs out of 41 fields could be 
associated with > 100 MeV photons. Thus, instead of 17 
LAT-detected GRBs over a two-year period, there might 
be 45 or 19% of all viewable GBM triggers with a lower 
limit of 25 or 11%. That is lower than any estimates 
prior to the Fermi launch. 

5. DISCUSSION 

By using optimal signal detection techniques, we have 
demonstrated that the existence of GRB high energy 
photon emission can be extended downwards in bright- 
ness by a significant factor. Increasing the dynamic range 
over which such phenomena can be measured puts a num- 
ber of constraints on physical models that can describe 
these events. As we have shown, these statistical meth- 
ods can set interesting limits on the behavior of ensem- 
bles of events at lower levels than is possible for single 
isolated bursts, raising serious questions about the ho- 
mogeneity of the GRB phenomenon. After another two 
years in orbit, the Fermi data will increase enough to al- 
low more careful re-examination of the initial conclusions 
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presented here. 




The most surprising aspect of this work is the very 
small number of GRBs that could be positively identified 
with high energy emission as a result of the substantially 
lower fluence thresholds. Almost any model would antic- 
ipate an increase in number proportional to the inverse 
of the ratio of detection thresholds raised to the power of 
unity or greater. Monte Carlo estimates indicate that our 
detection efficiency was of the order of 50% for the two 
events reported. This does not come close to explaining 
the apparent deficit of low-intensity events compared to 
the number of bursts similar to GRBs 080916C, 090510 
and 090902B . We look forward to pursuing this prob- 
lem during this precious and unique period when both 
Fermi and Swift are working magnificently. 



0.1 1.0 10.0 100.0 1000.0 
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Fig. 4. — Plot of high energy (LAT) fluences vs. low energy 
(GBM) fluences for the bursts considered in this paper. The six red 
points indicate the events previously reported by the LAT team. 
The two green points were identified by the matched filter tech- 
nique. The solid magenta line is a linear fit to the six LAT-detected 
bursts, 080916C, 090323, 090328A, 090510, 090902B and 091003. 
The dashed lines above and below indicate LAT fluences ten times 
greater and smaller. The blue point is the co-added sum of nine 
GRBs with GBM fluences comparable or greater to GRB 091208B, 
ie. greater than 5 /^ergs/cm 2 . The abscissa is the sum of all 9 GBM 
fluences and the ordinate is estimated from the total of all event 
class 2 & 3 photons in the corresponding burst fields. This indi- 
cates a significant population with high energy fluences below the 
approximately 1:1 fluence ratio represented by the diagonal solid 
magenta line. 
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a Swift trigger number if numeric; otherwise satellite name. 
b GRB Right Ascension. 
c GRB Declination. 

^ Angle between local orbit zenith direction and GRB. 
c Borcsight angle of GRB with respect to the LAT principal axis. 
GRB Galactic latitude. 

s Total number of LAT-detected photons over 250-s interval excluding those beyond 105° zenith cut. 
h 711, "2 and 713 arc the number of photons within 10.5° of GRB direction in 47.5-s window. 
1 Total event matched filter weight including the sharing parameter, £ . 

J GBM fluence, Sqbm, calculated as described in the main text. * indicates an entry with no available GBM fluence measurement; the value is estimated from the BAT fiuence by 
multiplying by a factor of 1.67 as derived from a joint sample of BAT and GBM measurements. Most other fluence values were taken from Table 2 in Guetta et al. 2010. The exceptions 
were obtained from the following references: GRB 080905A, Bissaldi et al., 2008; GRB 081118, Bhat ct al., 2008; GRB 090510, Guiriec et al., 2009; GRB 090518, Paciesas, 2009; GRB 
100206A, von Kicnlin 2010. 

k LAT fluence, Slat, calculated as described in the main text. 



